Min-protein oscillations in Escherichia coli with spontaneous formation of 
two-stranded filaments in a three-dimensional stochastic reaction-diffusion model 
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We introduce a three-dimensional stochastic reaction-diffusion model to describe MinD /MinE dy- 
namical structures in Escherichia coli. This model spontaneously generates pole-to-pole oscillations 
of the membrane-associated MinD proteins, MinE ring, as well as filaments of the membrane- 
associated MinD proteins. Experimental data suggest MinD filaments are two-stranded. In order to 
model them we assume that each membrane-associated MinD protein can form up to three bonds 
with adjacent membrane associated MinD molecules and that MinE induced hydrolysis strongly 
depends on the number of bonds MinD has established. 

PACS numbers: 87.17.Ee, 87.16.Ac, 87.16.Ka 



A division site in the rod-shaped bacterium E. coli is 
determined by the location of the FtsZ-ring 1] . Two ma- 
jor factors known to be important for placement of the 
FtsZ ring are nucleoid occlusion and Min-protein oscil- 
lations Q. Nucleoid occlusion restricts possible division 
sites to regions void of DNA — near the center and poles 
of the cell — while rapid pole-to-pole Min oscillations 
exclude poles as the possible division site 0, 0] . 

The Min system consists of three proteins: MinC, 
MinD, and MinE. MinD and MinE proteins generate 
pole-to-pole oscillations, while MinC proteins are being 
recruited to the membrane by MinD and hence follow 
the same oscillatory pattern Whereas MinC inhibits 
polymerization of FtsZ |(j , pole-to-pole oscillations pre- 
vent asymmetric cell division. MinD proteins in the 
ATP-bound form (MinD: ATP) attach to the membrane 
and presumably form two-stranded filaments arranged 
into a heli x El H MinE proteins function as ho- 
modimers [lOj; they attach to the membrane-associated 
MinD:ATP where they induce ATP hydrolysis, releasing 
subsequently MinD:ADP, MinE, and phosphate into the 
cytoplasm. The released MinD:ADP cannot bind to the 
membrane, until a nucleotide exchange takes place [j. 

There are several analytical 
and stochastic models [l8| that successfully reproduce 
Min oscillations. In this work we introduce the three- 
dimensional stochastic model. Our model, in contrast to 
others, takes into account a finite size of Min molecules 
and their spatial organization on the membrane. 



I. THE MODEL AND SIMULATIONS 

The shape of the bacterium E. coli is approximated by 
a cylinder of length H and radius R with two hemispheres 
of radius R at either end (two poles of the bacterium), 
giving the total length L = H + 2R (Fig. . Experimen- 
tally observed oscillations of MinC /MinD /MinE proteins 
between two poles are modeled using only MinD and 
MinE proteins. All interactions included in the model 
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FIG. 1: Geometric shape used to model the bacterium E. 
coli: cylinder of length H and radius R with two hemispheres 
of radius R at either end. The origin of the coordinate sys- 
tem (x, y, z) is placed in the center of the bacterium. Two 
distances from the membrane, d m i n and d ma x, are used to 
define the region near the membrane (d<dmin) and the re- 
gion far from the membrane (d > d max ), respectively. For 
details see text. All parameters shown in the figure are scaled 
equally, except the parameter dmin- 



take place simultaneously. However, for a particular 
MinD these interactions occur in four successive stages 
(Fig. |2l, which extend those proposed by Huang et al. 
|14j by taking into account the spatial organization of 
Min proteins on the membrane. 

(1.) Cytoplasmic MinD:ATP freely diffuses; When 
near the membrane it tries to attach to it in two ways: 
(i) either independently of other MinD: ATP molecules 
already attached, (ii) or it tries to became a part of a 
double chain of MinD:ATP molecules (two-stranded fila- 
ment) already formed on the membrane. 

(2.) MinE freely diffuses through cytoplasm. It does 
not attach to the membrane nor cytoplasmic MinD. 
However, MinE can attach to the membrane-associated 
MinD: ATP forming a MinE-MinD:ATP complex. 

(3.) MinE protein in the membrane-associated MinE- 
MinD:ATP complex stimulates detachment of the com- 
plex from the membrane by inducing ATP hydrolysis, 
releasing subsequently MinD:ADP, MinE, and phosphate 
into the cytoplasm. 

(4.) The MinD:ADP complex cannot attach to the 
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FIG. 2: Schematic representation of four stages in 
MinD/MinE proteins dynamics. (1.) MinD:ATP binds to the 
inner layer of the cytoplasmic membrane; (2.) MinE binds 
to the membrane-associated MinD:ATP; (3.) MinE induces 
ATP hydrolysis; releasing subsequently MinD:ADP, MinE, 
and phosphate into the cytoplasm, (4.) MinD:ADP is con- 
verted back into MinD:ATP by nucleotide exchange. 

membrane until it is transformed back into MinD:ATP 
by nucleotide exchange. 

In our model each molecule is treated as a point- 
like particle, except when on the membrane where the 
molecule size is taken into account. 

The diffusion process in three-dimensional space is 
described using Smoluchowski dynamics 0, A 
molecule starts from a well-defined position r at time 
t and diffuses during a time At. Probability density for 
finding a molecule at time t + At at a position r + Ar is 
described by 

P (f+Ar,t + At) = G s (Ax)G s (Ay)G s (Az), (1) 



s = V2DAt, (3) 

where G s (Ax) is a normalized Gaussian distribution with 
deviation s (diffusion length), and diffusion coefficient D. 

The position of the molecule at time t + At is obtained 
by adding random displacement to the current position 
where distribution of random displacements obeys {Tjl. 
Generally, time steps At do not have to be kept con- 
stant. In our simulation we use adaptive At in order to 
focus computational effort on important time segments. 



Because diffusion and unimolecular reactions are the only 
processes that take place in the region far away from the 
membrane, one can use longer time steps in that region 
than in the region near the membrane where, in addi- 
tion, bimolecular processes occur. These two regions are 
defined using two free model parameters d m in and d m ax 
— characteristic distances from the membrane (Fig. [JJ. 
In the region far away from the membrane (d > d max ) 
time step used is significantly longer than time step used 
in the region near the membrane (d < d m i n ). In the 
transitional region (d min <d<d max ) time step is gradu- 
ally decreased when approaching the membrane, to avoid 
that molecules entering the region near the membrane 
diffuse too far, avoiding on their path bimolecular reac- 
tions. For the same reason, time step At in the region 
near the membrane has to be chosen such that condition 
s = \/2DAt <§; d m i n is satisfied. 

In our model, parameter d m i n is also used as the re- 
action radius parameter for all bimolecular reactions. 
Hence, the cytoplasmic MinD:ATP molecule can attach 
to the membrane only when it is in the region near the 
membrane. Probability for this reaction is given by the 
simple intuitive formula: 

At 

Pd = o D - . (4) 

Grain 

The probability is proportional to time step At — the 
longer you wait it is more probable for a reaction to take 
place — and inversely proportional to d m i n to ensure 
that the number of reactions taking place depends only 
on the reaction rate parameter <jd and not on the model 
parameter d m i n used to define the near membrane re- 
gion. If the reaction occurs, the molecule attaches to the 
membrane with random orientation. However, our model 
forbids this reaction to take place if the position where 
the molecule should bind is already occupied by another 
MinD:ATP. 

Additionally, cytoplasmic MinD:ATP can react with 
MinD ATP molecules already attached to the membrane. 
Experimental data suggest that MinD:ATP attached to 
the membrane polymerizes into two-stranded filaments 
0, B 13 • Lacking experimental data on the interaction 
between membrane-associated MinD:ATP molecules, we 
utilize the simplest assumption in which each MinD:ATP 
molecule can form up to three bonds with adjacent 
MinD ATP molecules (Fig. [2J). The probability for cy- 
toplasmic MinD:ATP to occupy any free attachment site 
that is within reaction radius (r < d m i n ) depends on the 
reaction rate and- 

At 2tt ,o , . 

PDd = o Di — \ V = —d min . (5) 

An attachment of MinE to the membrane-associated 
MinD: ATP complex can take place if molecules are 
within the interaction radius (r < d m i n ) and there is 
no MinE molecule already attached. The probability for 
this reaction is 

At 2 it - 

PE = (?E-y- ; V=—d mm . (6) 
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FIG. 3: (a) Orthogonal projection of the membrane- 
associated MinD proteins onto a plane parallel to the line 
connecting two poles of the bacterium. Each projected pro- 
tein is represented with a dot. Five time frames are shown. 
They refer to times 0, \T, §T, §T, T; where T » 80 sec is 
the period of oscillation obtained with parameters specified 
in the text, (b) A portion of (a) is enlarged to clearly show 
two-stranded filaments. 
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MinE protein in the membrane-associated MinE- 
MinD:ATP complex stimulates detachment of the com- 
plex from the membrane by inducing ATP hydrolysis. 
The probability for this reaction might depend on the 
number of bonds particular MinD:ATP has formed with 
its MinD:ATP neighbors, and we assume that the number 
of bonds established decreases the reaction probability. 

(i) 

Let a de ', i = 0,1,2,3 stand for the detachment reaction 
rate when the MinD molecule has i bonds established. 
Hence, 
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and the probability for this reaction is 
1 _ exp (-<7«At) 



W 

Pde 



(7) 



(8) 



The transformation of MinD:ADP into MinD: ATP by 
nucleotide exchange is treated as unimolecular reaction 
with reaction rate (J A ) DP ^ ATP Hence, the probability 
for this reaction during time interval At is 



ADP~>ATP 

Pd 



1 — exp i 



ADP^ATP 



At) 



(9) 



II. RESULTS AND DISCUSSION 

In our numerical simulations we have fixed parameters 
related to the geometry of the cell (Fig.QJ to i? = 0.5 /xm 
(one- half cell width) and L = 4 /im (cell length) . Diffu- 
sion constants for cytoplasmic MinD and MinE proteins 



FIG. 4: Histogram for the number of the membrane- 
associated MinD (MinE) proteins versus x-coordinate, for 
the same time frames as in Fig. [3] Size of the bin used is 
0.08 rim. Time average of the number of the membrane- 
associated MinD (MinE) proteins over three periods is given 
in the last row. 



that are used, are in good agreement with measured val- 
ues for E. coli proteins of similar size |2l| : 



D, 



De = 2.5 /im 2 /sec. 



(10) 



Reaction rate parameters are chosen such that experi- 
mentally observed oscillations are reproduced: 

(j id = 0.01 /im/sec, arid = 0.003 ^m 3 /sec, 
(7 E =0.02/xm 3 /sec, ^p^atp 



1/sec. (11) 



These parameters are similar to the parameters that 
Huang et al. 0] used in their analytical model. How- 
ever, they use only one hydrolysis rate parameter, while 



in our model there are four 
obey Q with the ratio: 
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(12) 



Other ratios have been tested also. It is found that it 
is essential to take significantly smaller than 
in order to generate oscillations, whose period primarily 
(and strongly) depends on the parameter a^Pj . In the case 

when all aty were taken to be identical, the oscillations 
could not be produced even if other parameters of the 
model were varied substantially. 
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FIG. 5: Dependence of the oscillation period on the total 
number of MinD and MinE molecules, and the cell length. 
Points marked by □ are obtained with parameters used to 
generate Fig. [3] In (a) and (b) parameters are varied one at a 
time while keeping the others fixed. In (c), instead of keeping 
the number of MinD and MinE fixed, their concentrations are 
fixed at values used in Fig. [3]while cell length is varied. 



In our simulation we use 4000 MinD molecules and 
1400 MinE homodimers, reflecting the in vivo situa- 
tion |22j. The two-stranded filament width and the 
MinD monomer length are fixed to 6 and 5 nm, respec- 
tively §]. The model is evolved in time with time step 
At = 4 • 10~ 5 sec for processes far from the membrane 
(dmax = 0.1 /im) and Ai = 8 • 10~ 7 sec for processes near 
the membrane (d m i„ = 0.01 /im). In the transitional re- 
gion time steps are gradually decreased when approach- 
ing the membrane. We have tested the simulation by 
significantly varying parameters d m i n and At and the 
same results were obtained. 

With these parameters we have reproduced pole-to- 
pole MinD/MinE oscillations (Figs.Eland0} with the pe- 
riod T «80 sec, which is compatible with experimentally 
observed range (30 - 120 sec) 0. Initially, all MinE and 
MinD arc placed in the center of the bacterium. Other 
initial distributions were tried (e.g., uniform distribu- 
tion), and the same type of oscillations always appeared 
after the transient period lasting up to one oscillation 
cycle. 

Distributions of the membrane-associated MinD/MinE 
proteins do not oscillate in phase — MinE distribution 
lags after MinD distribution (Fig. |3} • This phenomenon 
has been seen in experiments, and it was described as an 
oscillating MinE ring [2^. However, when time-averaged 
both distributions have a minimum in the middle of the 
cell (last row in the Fig.QJ which reflects distribution nec- 
essary for proper cell division. This minimum is experi- 
mentally observed only in the case of the MinD protein 
oscillation (ltjj. For the time-averaged MinE distribu- 
tion there are only model predictions and they disagree 
on this point; e.g., there are models which predict, in 
contrast to our model prediction, that the time-averaged 
MinE distribution has a maximum in the middle of the 

ceii [HEi. 

To confirm its robustness, the model was additionally 
tested for a variety of experimentally observed phenom- 
ena. Fu et al. [23J have found experimentally that the 
oscillation period increases with the cell length. The 
overexpression experiments reveal that the oscillation pe- 
riod increases with the amount of MinD, and decreases 



o 
o 

CO 



o 
o 

<M 



::: 



§M||p|i 

mm 



- .. R.-J ; ;.r ;■• ■ a ■ ..... 



-6 -4 -2 2 4 6 
x [|im] 

FIG. 6: Space-time plot of the membrane-associated MinD 
proteins in the filamentous cell (L — 15 /im ). Each MinD 
is represented with a dot. For the sake of clearness, we have 
reduced the number of dots by a factor of 100. Concentrations 
of MinD and MinE, plus all the remaining model parameters 
are fixed at values used to obtain results in Fig. |3] 
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FIG. 7: Membrane-associated MinD proteins oscillations; 
only a single time frame is given. Parameter values are the 
same as those used in Fig. [3J except and = 0.02 /im 3 /sec, 
o\d = 0.0005 /im/sec, and er£P = 11.0 sec^ 1 [parameters er£,', 

a de 'j °"de' were modified to obey (1121 1. MinD filaments are 
longer then those obtained in Fig. |3] 



with the amount of MinE All these phenomena 

are reproduced with our model (Fig. |SJ. Consistent 
with experiments in the case of filamentous cells 
the zebra-striped oscillation pattern is obtained sponta- 
neously (Fig. El, starting with uniform distributions of 
MinD and MinE. 

In our model, MinD proteins attached to the 
membrane predominately form two-stranded filaments 
[Fig.Efb)]. This is achieved by imposing the ratio |Q to 

(i) 

the parameters tr^ e which are responsible for dynamics 
of both formation and decomposition of two-stranded fil- 
aments. The imposed ratio strongly favors two-stranded 
configurations over a group of single molecules — the 
probability for a group of single molecules to be detached 
from the membrane is considerably greater than that for 
the same group of molecules, but in the form of the two- 
stranded filament. 

The filament appears as an alive object. It is degraded 
and rebuilt constantly. When it grows in size the build- 
ing process dominates over the degrading process. Both 
processes preferentially take place at the filament's end. 
MinD molecules located at the end of the filament can 
form one or two bonds with its neighbors, while other 
MinD molecules have probably established three bonds. 
Because of I112H it is more probable for MinE to detach 
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MinD molecules locate the filament's end. If the building 
process dominates over the degrading process, detached 
MinD molecules will probably be replaced with cytoplas- 
mic MinD: ATP molecules. 

However, as the concentration of MinD molecules at- 
tached to the membrane reaches its peak, the concen- 
tration of cytoplasmic MinD:ATP goes to its minimum. 
At that time the degrading process dominates over the 
building process. Cytoplasmic MinE molecules continue 
to attach to the MinD molecules of the two-stranded fil- 
ament. MinD released into cytoplasm by MinE is in the 
form of the MinD:ADP complex and cannot bind to the 
membrane. However, MinE released into the cytoplasm 
by the same process attach to free attachment sites on 
the filament, thus speeding up its decomposition. 

The average length of filaments obtained with our 
model depends on the model parameters, particularly on 
(JDd and (7d. If we increase the parameter aod and/or 
decrease the parameter , the probability for attaching 
MinD to the filament already formed on the membrane 



(regulated by (JDd) will increase with respect to the prob- 
ability for starting a new filament formation (regulated 
by (in)- Hence, the average length of filaments is in- 
creased (Fig 01. In order to keep the period of oscillation 
similar to the period for the case shown in Fig. |3| hydrol- 
ysis rate parameters were modified: a^J = 11.0 sec -1 , 
while the same ratio (1121 was obeyed. 

In conclusion, we have introduced 3D stochastic 
reaction-diffusion model to describe MinD /MinE dynam- 
ical structures in E. coli. In particular, our model spon- 
taneously generates two-stranded filaments using a few 
simple physical assumptions. 
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